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Abstract 

Dielectric nano-swithes made of the materials that exhibit piezoelectric and /or flexoelectric prop- 
erties with significant electro-mechanical coupling are considered. In this case, a nonuniform strain 
field may locally break inversion symmetry and induce polarization even in nonpiezoelectrics. At 
reducing dimensions to the nanoscale, the flexoelectric effect demonstrates the nonlocality of the 
dielectric materials and plays more significant role than piezoelectric effect. The flexoelectric ef- 
fect is included into consideration via additional term coupling strain gradient and polarization 
in the electric enthalpy density. The equations of motion of the improved Euler-Bernoulli and 
Timoshenko beam models, and 2-D plate theory have been obtained. 



Electronic address: marina. toropova@utoronto.ca 



I. INTRODUCTION 



Electromechanical nano-switches as integral parts of nanoelectromechanical systems 
(NEMS) find numerous technological applications as, e.g., mass memory storage, high- 
frequency electrical switches, and mass or force sensors. Such high technological applications 
demand combined efforts of engineering and science through modelling and simulations. 
Hence, to design state-of-the-art nanotechnological devices with predetermined character- 
istics, engineers need in addition to experimental techniques not only new formulas and 
computational methods but also improved electromechanical models. The new models must 
account for the nonlocal properties of the materials and new physical phenomena. 

The nonlocality appears due to the noticeable role of interatomic forces in nano-objects. 
In nonlocal theories, the constitutive equations take into account microstructure of real 
materials and microscopic interaction length between, e.g., molecules in a lattice. From the 
physical point of view, the nonlocality exhibits through special size- dependent effects such as, 
for example, flexoelectricity - induced polarization due to the strain gradient. It is well known 
that conventional continuum theories are size-independent and therefore cannot be applied 
automatically for the analysis of NEMS devices. The most suitable tools for their analysis 
are atomic and molecular models, but they are restricted by their computational capacity. 
One of the ways to resolve this contradiction consists in the employment of improved, size- 
dependent classical theories [1]. 

In the present work, we consider both piezoelectric and flexoelectric cases. In the latter 
case, a nonuniform strain field may locally break inversion symmetry and induce polarization. 
To take into account this effect, similar to [1], we assume that the electric enthalpy density 
depends not only on strain, electric field, and polarization but also on strain gradient. 
Then, we replace the polarization by its linear representations through the electric field and 
strain. Using Hamilton's principle, we obtain the equations of motion for Euler-Bernoulli 
and Timoshenko beams as well as for 2-D plate. The higher order terms in the strain 
from the electric enthalpy density give additional contribution to the bending rigidity of 
the beam compared to classical solution. The analysis shows that taking into account the 
strain gradient increases the elastic characteristics of the nano-switch considered in this 
work. From the obtained formulas it is seen that at reducing dimensions to the nanoscale 
the flexoelectric effect plays more significant role than piezoelectric effect. In the dynamic 
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case, the Euler-Bernoulli beam model provides overestimated frequencies. Two-dimensional 
model allows us to evaluate the electric potential accumulated in the nano-switch due to 
induced nonuniform strain. 



II. GOVERNING EQUATIONS 



Electromechanical nano-switch may be simulated as a dielectric cantilever nanobeam 
(Fig. 1). The material of the beam may have tetragonal or cubic symmetry or in another 
words it may be piezo- or nonpiezoelectric. 

n n 

First, let us consider flexoelectric properties of the beam. In [1|]-|4| the fourth-order 
flexoelectric tensor is introduced in two different ways: as fiijkiUj^i 



Pi — ZoVlmEm + ZljkSjk + f^ijklUiJk 

r~| 

or as fijki in the additional term in the thermodynamic potential |4|]: 

fijkl-Pi^jjkl- 



(1) 



(2) 



Here, SV,-, ey&, rji m , are the components of the strain, dielectric, and relative permittivity 
tensors, respectively, Pi, Ui, E m , are the components of polarization, displacement, and 
electric field vectors, respectively, e a is the permittivity of vacuum. 

In jlj] there is a structure of the flexoelectric coefficients for a cubic material (in matrix 
notation): 
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According to |l[ and [6j, the electric enthalpy may be presented in the form 



H(Sij) Pi, Pi t ki U 



k, Uijk 



-e (p,i(p,i + tp,iPi, 



(3) 



— ClijPiPj + -^bijklPj,iPl,k + ^CijklSijSkl 



W(Sij, Pi, Pi t k, Ui,jk) 
"^dijkSjkPi -f- &ijklPj,iSkl fijkl-Pi^'jykl 9ijk-Pi-Pk,j j 



(4) 




FIG. 1: Cantilever beam(initial configuration). 



where W is an energy density of deformation and polarization, ay, b^ki, c^z, d^ e^h 
fijki, Qijk are the components of reciprocal dielectric susceptibility, polarization gradient- 
polarization gradient coupling tensor, elastic tensor, piezoelectric tensor, polarization 
gradient-strain coupling tensor, flexoelectric tensor, and polarization-polarization gradient 
coupling tensor, respectively, if is the potential of electric field. 

Since the components of the stress tensor and the electric field vector Ei = —<p i are 
expressed as 

dW dW 

tj ~ ds~/ 1 ~ ~ws 

we have 

T{j = CijklSkl + dijkPk + CijklPk,l 
Ei = CLijPj + dijkSjk + fijklUijk (5) 

Compare (JTj), and the second equation of (|3j), we conclude that flexoelectric tensor faki 
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has 21 nonzero components: 
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In this work, we do not consider the terms with bijki, tiju, and coefficients although it 
might be done without difficulties. 

Using the equations (JS])-© we derive the equations of motion for a dielectric nanoswitch. 

Here we consider a cantilever switch which may be presented as a Euler- Bernoulli or 
Timoshenko beam or 2-D narrow plate (Fig. [TJ). The load applied to the free edge of the 
cantilever simulates the case of nonuniform strain. 

III. EULER-BERNOULLI BEAM MODEL 

Within the framework of the Euler-Bernoulli model, the displacement can be written as 

ui(x, y,z,t) = u ~ — zw (x, t), ii2(x, y, z, t) « 0, 

u 3 (x,y,z,t) » w(x,t), (6) 

Hence, based on the equations (EJ), the strains are equal to 

Sn = e x = -zw", S 2 2 = £ y = 0, S 33 =e z = 0, 
S12 = S 2 i = ixy = S 13 = S 31 = j xz = S 23 = S 32 = jyz = 0. (7) 

We assume that the electric field acts only in z— direction. Then, from (jSJ) we can express 
E z as 

— E 3 = —E z = d 33 P 3 + d 3 iE x + /1133W + fl313U zx + f3113W,xx = 

033^2 - (d 31 z - f 113a )w" 
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or 

P 3 = P Z = —(-E g + (/nss + d 31 z)w"), (8) 

where the low letter indexes after a comma mean the derivation with respect to corresponding 
coordinates. 

After substituting the equations (|7|), (|S]), into (j3J) and (jlj) and ignoring the polarization 
gradient, we have 

H = \ c ll z2 { w " f + 7T—(Vz + (Al33 + d 31 z)w") 2 - 
I A(X 33 

// 

d 31 w" (<p' z + (/ 1133 + d 31 z)w") - /n33 (cp' z + (/H33 + rf 3 iz)w") 

033 033 

-V^ 2 -— (^ + (Al33 + 4lV) (9) 
^ «33 



5 [\T k -H + W d )dt = 



The equation of motion of the beam are derived via Hamilton's principle 

-t 2 



(10) 



where 5(-) denotes the first variation, Tjt is the kinetic energy, and Wd is the work done by 
the external forces and moments. 

The kinetic energy of the beam by using the equation (J6]) can be written as 

Tk = \ I Q{u} T {u}dV b = \l q[w 2 + z 2 (w) 2 ]dV b = 
1 Jv b 1 Jv h 

\ e [ [(w) 2 A + I(w') 2 }dx, (11) 
1 Jo 

where I = J A z 2 dA, A is the area of the beam's cross section, V& is the volume of the beam, 
q is the mass density of the nano-switch's material, and the length of the beam is equal to I. 
If the beam is under the transverse force q, the work done is equal to 

W d = [ qwdx. (12) 
Jo 

Performing variation in (fit)]) and using flU]), (jTTj) . and (112"]) . we obtain the equations 
describing the electromechanical behavior of the beam 

d 2 T f 2 

^ (ciiJ _^_/ll33 A) _ 
a 33 a 33 

1 fP 1 f)^W 

ifnssQt 1 + d 3l Mf) + Q Aw + Ql irTT -, = q, (13) 



a 33 dx 2KJLLM ^ z ox z ' e dx 2 dt 2 



where Qf and Mf are the electric transverse shear force and the electric bending moment, 
respectively: 

Qf= [ E z (x,z)dA 

J A 

Mf = [ zE z (x,z)dA 

with the associated boundary conditions 
at x = 



A 



w — w' — 0, (15) 



( Cll J f -h*A)w" - d -^M? - ^Qt l = 0, 

^33 a 33 a 33 a 33 

(cu/ _ §£ _ ^33^)^'" _ ^i|-M^ - (16) 

Gt33 «33 «33 OX 

/U33 d_ QeJ = ^ 



a 33 dx 



and at z = ±h/2 



<P, = 0, (17) 



where P* is a force applied to the right edge of the cantilever. 

The equations ( IT3l) . ( JT4l) and the boundary conditions ( 1T51) - ( IT71) couple the mechanical 
displacement w and the electric potential ip. 

From the equation f lT3|) it is seen that piezoelectric and flexoelectric effects increase the 
bending rigidity of the beam (033 < 0) and for the beams with large cross-section area 
and consequently large its moment of inertia, the term Af 2 133 /a 33 plays insignificant role. 
However, for small cross-sectional dimensions, when 

f 2 d 2 h 2 

J 1133 ^ („ U 31 \ u 

> [Cn Jttt, 

033 «33 tl 

where h is the thickness of the beam, the flexoelectric effect becomes noticeable. 
The normalized Young modulus is 



Y 



^31 12/ 2 133 



C11CI33 cna 33 /i 2 



It is seen that at h — > the flexoelectric term plays the dominant role in the bending rigidity. 

The classical formula for the maximum delection of a cantilever beam w max = PJ 3 /3Y 
may be used to find the coefficient /1133 from experimental results: 

_ / 1 ' h? pj~ } 

J1133 — \ C HT7T — 7T7T 

y a 33 1 2 3bhw max . 
It is obvious that the dispersive relation is 

u 2 g(A - Ik 2 ) - k\ Cll I -^L- G™A) = (18) 

^33 a 33 

where u is a frequency and k is a wavenumber. 



IV. TIMOSHENKO BEAM MODEL 

In this model we assume that the displacements are presented as 
m(x, y, z, t) « -zij)(x, t), u 2 (x, y, z, t) « 0, 

u 3 (x,y,z,t) « w(x,t), (19) 

where ty is the transverse displacement of the points of the centroidal axis(y = z — 0), and 
^ is the rotation of the beam cross-section about the positive y— axis. 
From f|T9|) . the strains are equal to 

e x = -zif)', e y = 0, e z = 0, 

Ixy = lyz = 0, J xz = W - if). (20) 

Since the narrow beam deflects in x — z plane, we consider only two components of the 
electric field E x and E z . From (E]) they can be expressed as 

—E x = a\\P x + di^ xz + fmiu,xx 

—E z = a^P z + d,3iE x + fll33U jXZ + fl313U tZX + f3113W tXX: 



or 



-.Ej = oiiP x + d 15 (w -ip)- fuuzip' 

-E z = a 33 P z - d 31 zip' - /ii 33 ^' + fiais{w" - if)'), 
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and hence 

P x = —{-E x - d 15 (w' + fimzip"), 
an 

P Z = —{~E Z + d 31 ZlJ>' + /H33V>' - W™" - V')), (21) 
^33 

After substituting the equations (1201) . (12T]) into ([3]) and (j4]) and ignoring the polarization 
gradient, we have 

1111 

H = -a u P 2 + -a^P 2 z + -ch^ + t^sstL + di 5 ^ xz P x 

+d 31 e x P z + fimPiu xx + P 3 (/ii3 3 M^ + fnnUzx + hmWxx) - 
l^oiv'x) 2 - \^(v' z ? + v' x Px + ¥>' z Pz = 

1 '(fx - d ^( w ' - + /llll^") 2 + 



2«n 

^— + d 3 i^' + /ii33^' - /i3i 3 (w" - V')) 2 + 

1 / 1 / i\i i a 
-diz 2 ^') 2 + -c 55 (w - ip) 2 + d l5 {w' - ip) — (<p' x - d l5 (w' + /mizV") 

2 2 an 

d 31 zip , + + j U33 ^' _ f 1313 ( w " - ijj')) - 

a 33 

hm- — (jp' x - di 5 (w' -i>) + /imV') 

G33 

/ii33 — {fz + d 3 i2V' + /nraV' - /l313(w" - -0')) + 

^33 

/l313 (V 3 ! + d 3 lZ1p' + fll 33 1p' ~ fl313{w" ~ 1p')) 

«33 

1 1 ' 

-^o(p'x) 2 ~ o^o(^l) 2 + — (<f' x ~ di 5 (w' - t/j) + fmiz^") + 

2 2 an 

/ 

— (<p' z + d 31 ztp' + f 1133 i/)' - fi 3 i 3 (w" - ip')) 
a-33 

The kinetic energy is equal to 



T k = I [ g{u} T {u}dV b = \ I g[zr + w 2 }dV b = \q / [/(^ 2 ) + A(w) 2 }dx, 
1 ->v b 1 Jv b 1 



The work done is described by the equation (fT2|) . 

Applying the Hamilton's principle similar to previous section, we obtain the equations of 
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motion for the Timoshenko beam 



kA(c 55 -^)(w"-^) + ^Ql l + 
a n an Ox 

— i-f^Q? ~ M/llSS^' " /iSlS^ - /'))) - ? = ^ 
033 ox 

(cnl -^1- { f 1133 + f 1313f kA)^" + kA(c 55 - d ^){w - 

CI33 CI33 flu 

an ox z an 

/ 3 . ,3 . ,di 5 d3i. / di 5 » 

£o) + ^( £ o) + ( 1 )^ w =0 

an 033 o-n 033 an 



with boundary conditions: 
at x = 



w = w — ip — tp — (f — 0, 



and at rr = /: 



k^(cte-^)(«/-^) + — <# + 

an an 
Ag^ _ kA{fu3 ^ _ Mw »> _ /))} 
033 0^ 



(c n 7 - - (/ll33 + /l313)2 fcA)^ - —Mf + 

0-33 O33 CI33 

/l!33 + /l313 , , , „ " —ix . /llll <9 d /• 



a 33 an dx a u 

kA[(f 1133 + / 1313 )V' - /W] + Of = 0, 

/1111/V" + = 0, 



y? = 0, 



and at 2; = ±h/2 



, 3 \ ' . /1133 

3 3 «33 



1313 



1313 



= 0. 



033 
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Here k = 5/6 is the shear correction factor and 



E x (x, z)dA. 



Ml 



Now, let us introduce new quantities 



zE x (x, z)dA. 



Cl 



\ 



„ <*li (/1133+/1313) 2 

Cl 1 2 

«33 £133^0 

Q 



r 



Then, the dispersion relation can be written as 



oo 4 - c/p§ + (c 2 s + ci)k 2 } + cicik* = 



or 



i((cf + ft* + 4) ± J-Mci + + %Y - cjcft' 

Z r V 



(22) 



The equation fl22|) presents the frequencies for two wave modes. The lower frequency relates 
to flexural wave and the higher frequency relates to shear wave. 

From the formulas ( )22|) and ( fTHI) it is seen that the eigenfrequencies increase if we take 
into account the piezoelectric and/or flexoelectric effects. The results of the comparison 
between the eigenfrequencies calculated within the framework of the Euler-Bernoulli and the 
Timoshenko beam models coincide with the conclusions related to single-walled nanotubes 
presented in 



V. TWO-DIMENSIONAL MODEL 

In this model, we consider two components of the displacement u — U\ and w = M3 and 
only three components of strain tensor acting in x and z direction 
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The electric field may be written as 



-E x = a n P x + dy^Tixz + fuuu xx + fi33i(u xx + w xz ) + /3311W31 
-E z = a 33 P z + d 31 u x + f U33 u xz + fm 3 {u 31 + w u ) + f 3333 w zz . 



Now, the polarization may be expressed as 



Px = — —{ l P x ~ d 15 {u z + w x ) - fu U u xx - f 133 i{u xx + w xz ) 
an 

/3311W31) 

Pz = — -(<f' z - d 31 u x - fi 133 u xz - fi 3 i 3 (u 31 + w n ) - f 3333 w Z: 
a 3 3 



In this case, the equation (j3J) takes the form 



H = ^a u P 2 + ~ a ™P z + \c xl e 2 x + \c 33 e\ + ~c 55 7L + c 13 e x e y + 

dl5lxzPx + d 31 e x P z + Pl{fnnU xx + Fi 331 (u zz + W xz ) + f 33 nW xz ) + 



^(/ll33«ra + fl313{u xz + W xx + f3333W zz ) ~ 
1 / 1 / / / 

2 £ o(y?L) 2 - 2 £ o(^l) 2 + Vx p x + f z p z = 

-—((f' x - d 15 (u z + W x ) - flUlU xx - fl 33 l{u zz + W xz ) - f 3 31lW xz ) 2 

2a n 

+ ~ {<P Z - d 31 U x - fu33U xz - f V3 l 3 {u zz + W xx ) ~ f3333W zz ) 2 + 

*Q>33 

111 

-c n ul + -c 33 w 2 z + -c 55 (u z + w x ) 2 + c 13 u x w z + 

— (lf' x - d 15 (u z + W x ) - flinUxx - fl33l(Uzz + W xz ) - f 33 uW xz ) ■ 

an 

(f'x + dis( u z + w x ) + f uu u xx + fi 3 3i(u zz + w xz ) + f 33 uw xz ) + 
1 



2a 3 3 



(<P Z ~ d 3 iU x - fl!33U xz - fl3is(u zz + W xx ) - f 3333 W zz ) 



I 1 / 1 / 

(<f' z + d 31 U x + fu 3 3U xz + /l313(w** + Wxx) + f3333^ zz ) ~ T^ofax) 2 ~ ly^iVz)' 



The kinetic energy is 



The work done is 



T k = \ I e [v? + w 2 ]dV b 

Wd = / (pu + qw)dx, 
Jo 
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where p is the longitudinal load and q is the transverse load. Via the Hamilton's principle, 
we have the equations of motion 

A 2 a 2 A 2 A A 

{ 31 \ / 15 \ / 15 \ / ^31 ^15 x 

u xx {c u ) + w 22 (c 55 ) + w xz {c 13 + c 55 ) - if xz { 1 ) + 

Cf33 a ll a ll O33 On 

/ dm r d 3 i /1331 d 3 i di 5 
w xxx \ — jmi J1313J H «i5 M zzz — w xzz { — J3333 H J1331J 

Oil «33 a ll ^33 a ll 
[<fxxx + fllllU xxxx + fl33l(u xxzz + W xxzz ) + /33IIW xxxzj 

an 

— ■ \fxzz + fllll u xxzz + fl33l( u zzzz + w xzzz) + f331l w xzzz] 

an 

f^. 133 i/*131- 3 

Wxzz + f\\33 u xxzz + fl313( u xxzz + W xxxz ) + f3333W xzzz ] — p = QU, 

a 3 3 

d 2 d 2 

w xx (c 55 -) + W ZZ C 33 + U xz (c 13 + C55 -) (23) 

a ll a ll 

^15 ^15 r /dsi r di 5 

Pxx J33ll w xxz + U xxx { J1313 JllllJ — 

On On ^33 a ll 

[<fxxz + /ll33« xxxz + fl313{u xxxz + W 

xxxx ) + /3333W xxzz + 

033 

[fxxz + /llllW + fl33l(Uxzzz + W xxzz ) + /3311W xxzz I 

an 

^^-[fzzz + fll33Uxzzz + fl313(Uxzzz + W xxzz ) + / 3 333W 2222 ] - q = QW , 
033 

/ 3 . ,3 . ,4 . 4\ ^15 

<Pxx{ £0) + £0) - U xz { 1 ) W xx - 

an 0*33 a ll a 33 °11 

/llll _ /3333 _ / /l331 /l313 x / /l331 /3311 /l313 x _ n 

Uxxx W zzz Uxzz\ r ) ~r W xxz y I r ) U 

ail «33 a ll ^33 a ll a ll «33 

with the boundary conditions presented for the brevity sake in the form of partial derivatives 
of the electric enthalpy function as: at x = 

u = u' x = w = w' x = <jd = 0, (24) 



at x — I 



*■ 



dH _ d_^dH_^ _ _ p 

du x dx du xx 2 dz du xz 

ou xx 

^IL _ ^_(HL) _ I _^__^L o (25) 

dw x dx dw xx 2 dz dw xz 

dw xx 
dH -0 
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and at z — ±h/2 

dH d ,dH , 1 d dH 



du z dz du zz 2 dx du 3 
dH =0, 



0. 



du z ,. 

dH d , dH 4 I d dH 



dw z dz dw zz 2 dx dw xz 
dH 



(26) 



dw zz 
dH 



= 0, 
0. 



dip z 

If the flexoelectric coefficients are equal to zero, the equations (|23|) coincide with the equa- 
tions for plates presented in 

Numerical solution to the system (I23p with the boundary conditions §M§ - ( 12^|) allow us 
to find the deflection of the nano-switch, the strains, the stresses, and the electric potential 
accumulated in the cantilever. 



VI. CONCLUSIONS 



In the present work, the nonlocal properties of the dielectric materials at nanoscale have 
been taken into account through the flexoelectric effect. Based on the Hamilton's principle 
and the electric enthalpy density with additional term describing the coupling between 
the polarization and the strain gradient, the equations of motion for Euler-Bernoulli and 
Timoshenko beam models as well as 2-D plate model have been derived. These equations, 
may be used to analyse the static and dynamic behavior of cantilever nano-switch. The 
formula connected the flexoelectric coefficient and the maximum deflection of the cantilever 
has been presented. 

Two dimensional model allows us to analyse the electric potential accumulated in the 
nano-switch due to the flexoelectric effect. 
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